Zero-modes of the QED Neuberger Dirac operator* 
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We consider 4d compact lattice QED in the quenched approximation. First, we briefly summarize the spectrum 
of the staggered Dirac operator and its connection with random matrix theory. Afterwards we present results for 
the low-lying eigenmodes of the Neuberger overlap-Dirac operator. In the strong coupling phase we find exact 
zero-modes. Subsequently we discuss possibly related topological excitations of the U(l) lattice gauge theory. 



1. Spectrum of the QED Dirac operator 

The spectrum of the QCD Dirac operator 

iD + im = \ ■ in a chiral basis (1) 

\^ 1 ' mi J 

is related to universality classes of random matrix 
theory (RMT), i.e. determined by the global sym- 
metries of the QCD partition function, for a re- 
view see [0 . In RMT the matrix T in Eq. (|) IS re- 
placed by a random matrix with appropriate sym- 
metries, generating the chiral orthogonal (chOE), 
unitary (chUE), and symplectic (chSE) ensemble, 
respectively. For SU(2) and SU(3) gauge groups 
numerous results exist confirming the expected 
relations. 

We have investigated Ad U(l) gauge theory de- 
scribed by the action 5'{C/i} = X]p(l~cos0p) with 
Ui = Ux.j^ = e^p{idx,^) and Op = O^,,^ + B^+jx,^ - 
Ox+o,,, - dx,u iy^pi). At /3c « 1.01 U(l) gauge 
theory undergoes a phase transition between a 
confinement phase with mass gap and monopole 
excitations for (3 < (3c and the Coulomb phase for 
[3 > Pc- In the whole Coulomb phase the pho- 
ton is massless and thus there ought to exist a 
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continuum theory, asV^oo, everywhere. 

We were interested in the relationship between 
U(l) gauge theory and RMT across this phase 
transition. The Bohigas-Giannoni-Schmit conjec- 
ture j^] states that quantum systems whose clas- 
sical counterparts are chaotic have spectral fluc- 
tuation properties, measured, e.g. by the nearest- 
neighbor spacing distribution P{s) of unfolded 
eigenvalues, given by RMT, whereas systems 
whose classical counterparts are integrable obey 
a Poisson distribution, P(s) = exp(— s). With 
staggered fermions it was shown in Ref. ||^ that 
the Bohigas-Giannoni-Schmit conjecture holds 
for lattice QED both in the confined phase (as 
for the SU(2) and SU(3) gauge groups) and also 
in the Coulomb phase, whereas free fermions yield 
the Poisson distribution. In Ref. [Q this investi- 
gation was continued with a study of the distri- 
bution of small eigenvalues in the confined phase. 
Excellent agreement was found for the conjecture 
that the microscopic spectral density 

Psiz)^ lim —p(—\ (2) 

is given by the result for the chUE of RMT, 
which also generates the Leutwyler-Smilga sum 
rules [0 . Here S is the modulus of the chiral con- 
densate, which follows from the smallest non-zero 
eigenvalue via the Banks-Casher formula. The 
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Figure 1. Density p(A) of small eigenvalues for 
the S'^ X 6 lattice at /? = 0.9 (left plot) and at 
f3 — l.l (right plot). A non-zero chiral condensate 
is supported in the confinement phase. 

quasi zero-modes which are responsible for the 
chiral condensate build up when we cross from the 
Coulomb into the confined phase. For an 8'^ x 6 
lattice [|, Fig.! compares on identical scales den- 
sities of the small eigenvalues at (3 = 0.9 (left plot) 
and at /3 = 1.1 (right plot), averaged over 20 con- 
figurations. The quasi zero-modes in the left plot 
are responsible for the non-zero chiral condensate 
S > via the Banks-Casher formula, whereas no 
such quasi zero-modes are found in the Coulomb 
phase. For M SU(2) and SU(3) gauge theories 
a general interpretation is to link zero-modes to 
the existence of instantons. As there are no con- 
ventional instantons in 4d U(l) gauge theory (an 
analogous case exists in 3d QCD [^), we Q de- 
cided to study the physical origin of the U(l) 
quasi zero modes in more detail. 

2. Exact zero-modes 

Via the Atiyah-Singer index theorem the topo- 
logical charge is mapped on the number of 
fermionic zero- modes of the Dirac operator. It 
is thus desirable to use a Dirac operator which 
retains chiral symmetry and exhibits exact zero- 
modes. Therefore, we employed in Ref. the 
Neuberger overlap-Dirac operator ||^,^ 
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Here 75i/^(— m) is the usual Wilson-Dirac oper- 
ator on the lattice and e the sign function. The 
mass TO has to be chosen to be positive and well 




Figure 2. Exact RMT probability densities (|) 
for the lowest non-zero eigenvalues compared with 
histograms from our data. 



above the critical mass for Wilson fermions but 
below the mass where the doublers become light 
on the lattice. 

Our first question was whether the QED 
overlap-Dirac operator will exhibit exact zero- 
modes at all. To answer it, we have analyzed 
configurations on L'^ lattices at (H — 0.9 in the 
confined phase and at f3 — 1.1 in the Coulomb 
phase. With an overrelaxation/heatbath algo- 
rithm we generated 500 configurations per lat- 
tice of linear size L = 4, 6 and 8. After ther- 
malization, the configurations were separated by 
100 sweeps, with each sweep consisting of three 
overrelaxation updates of each link, followed by 
one heatbath update. On each configuration the 
lowest 12 eigenvalues of the overlap-Dirac opera- 
tor (||) were calculated using the Ritz functional 
algorithm with the optimal rational approxima- 
tion of Ref. @ to e {Hyj{m)) and to set to 2.3. 

In the confined phase exact zero-modes of the 
operator (||) were indeed found. The highest de- 
generacy observed was = 3. No zero-modes 
were found in the Coulomb phase. For each de- 
generacy V, chiral RMT predicts the distribution 
of the lowest non-zero eigenvalue Amin in terms of 
the rescaled variable 

^ = SV^A,„i„, (4) 

where V is the volume of the system and S the in- 
finite volume value of the chiral condensate (V'V') 
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Figure 3. Decomposition of the spectrum p{X) (normalized to the number of eigenvalues) of the staggered 
Dirac operator into a monopole part Pmon(A) and a photon part /Ophot(A), averaged over 500 configurations 
on a 4'* lattice. 



determined up to an overall wave function renor- 
malization. For the U(l) gauge group the unitary 
ensemble applies and the RMT predictions for the 
V = 0,1,2 probability densities Pi^j^n(z) are |l| 



(5) 



The chiral condensate is related to the expecta- 
tion value of the smallest eigenvalue. For degen- 
eracy v we have 



(6) 



where the result is supposed to be independent of 
v and 



(7) 



Using analytical and numerical integration the 
expectation values {z^'^'') are easily calculated. 
The values S'-'^-' follow then from Eq. (||) and, for 
the various lattice sizes and j^'s, the results which 
we obtain from our data are consistent with one 
another. 

Figure || displays the exact RMT probability 
densities pl^}^ of Eq. (||) and the corresponding 
histograms from our data. It shows that the his- 
tograms follow the shift of the RMT probability 
densities and their general shape. The high peak 
of the v = 2 histogram is interpreted as a statis- 
tical fluctuation due to the low statistics we have 
for this case. 



3. Topological objects 

The topological structure of the U(l) gauge 
theory is not as extensively studied as the non- 
Abelian case. Nevertheless, with torons [|lO| , 
monopole solutions and Dirac sheets |l2| a 
number of topological objects are known for the 
U(l) gauge theory. All these topological configu- 
rations have in common with the zero-mode de- 
generacy that they are turned on at the phase 
transition from order (Coulomb) to disorder (con- 
finement). However, in the confined phase we 
found no detailed correlation between any of the 
topological phenomena and the zero-mode degen- 
eracy of the overlap-Dirac operator. Related to 
this may be that we found the zero-mode suscep- 
tibility to decrease, possibly approaching zero for 

In contrast to these findings there is an inter- 
esting observation that ordered Dirac sheet con- 
figurations give rise to zero-modes of the overlap- 
Dirac operator, with the number of zero-modes 
being equal to the number of Dirac sheets |§,|l3| . 
One might understand this from the fact that a 
Dirac sheet is a 2c? gauge configuration that con- 
tains unit topological charge in the 2d sense, kept 
constant in the two orthogonal directions. 

To investigate possible correlations between the 
existence of zero-modes and topological objects 
further we have, following Refs. p4|-p^, factor- 
ized our gauge configurations into monopole and 
photon fields. The U(l) plaquette angles 9x,fj,ii are 
decomposed into the "physical" electromagnetic 
flux through the plaquette 9x^1, and a number 
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'nix,^v of Dirac strings passing through the pla- 
quette 

6x,tiy — Ox^fiu + 2n nix, f^i, , (8) 

where ^^^^jy G tt, +7''] and m^^f^^ 7^ is cahed a 
Dirac plaquette. Monopole and photon fields are 
then defined in the foUowing way 



G:r rr' 9',, 



(9) 



(10) 



Here d' acts on x', the quantities rrix^^i, and 
0x,tj.i' are defined in Eq. (||) and Gx,x' is the lat- 
tice Coulomb propagator. One can show that 
9x,tj. = 6*"'°" + 9P^°^ is up to a gauge transfor- 
mation identical with the original 6x,i_i defined by 
Ux,^ exp(i6'^,^). 

We found that both zero-modes and near zero- 
modes lie solely in the monopole part of the 
gauge field and are completely absent in the pho- 
tonic field. Using periodic boundary conditions 
in space and anti-periodic boundary conditions in 
time , this was seen both for the overlap-Dirac 
operator and for the quasi zero- modes of the stag- 
gered Dirac operator. As an example we show in 
Fig. 1^ the entire spectral density of the staggered 
Dirac operator in the original (full) gauge fields 
and in their monopole and photon part, respec- 
tively. 

At the moment we are accumulating statistics 
to perform a decomposition into different topo- 
logical sectors and to obtain an analogous analy- 
sis concerning topological objects as we did with 
the original U(l) field above It is of further 
interest to study space-time correlations between 
different topological objects. We pose the ques- 
tion of the existence of local correlations between 
the topological charge density and the monopole 
density. It is also desirable to calculate spatial 
correlations between the topological objects and 
the density ip'^ip{x) of the exact zero-modes of the 
overlap-Dirac operator. 

4. Cooling 

In non-Abelian theories it turned out that some 
cooling or smoothing procedure is unavoidable to 
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Figure 4. Topological charge of artificially con- 
structed configurations on an 8^ lattice. 



gain integer- valued topological charges. For our 
Abelian case we studied the value Q obtained by 
integration of the topological density 



Q 
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327r: 



(11) 



using the hypcrcube definition and plaquette def- 
inition on the lattice jl7|. In Ad U(l) theory with 
periodic boundary conditions one can construct 
an artificial configuration with integer topologi- 
cal charge Q Mm 



Uiix 

U2{X 
U2{X 
U3{X 

Ui{x 
Uiix 



exp(iu!iax2) 

1 {x2 = a,2a,...,{N - l)a) 

exp(— icjiaTVxi) {x2 — Na) 
exp (107202:4) 

1 {x4 = a,2a, ...,{N ~ l)a) 

ejip{~iuj20-N X3) {x4 — Na), (12) 



where N is the linear lattice extent, w^a^ = ^j^, 
rii = —2, —1, 0, 1, 2, and Q = nin2- These 
gauge fields consist of ni Dirac sheets in the 1- 
2 planes and n2 Dirac sheets in the 3-4 planes. 
Such configurations possess nearly integer-valued 
Q's on a finite lattice, as seen from Fig. ^, and the 
overlap-Dirac operator has Q exact zero-modes. 

We have first heated and then cooled those 
smooth configurations in order to inspect whether 
they return to their original value of Q. In Fig. o 
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Figure 5. Cooling history of the topological 
charge of the Q = 1 configuration after increasing 
number of updates a,t (3 — 0.9. Integer topological 
charge is lost after i 10 updates. 



we demonstrate the effect for the Q = 1 configu- 
ration which was updated several times with the 
Monte Carlo code at P = 0.9 before 400 cool- 
ing steps were applied. We employed a system- 
atic cooling procedure by U^(x) VU^{x), with 
V = and A = U,ix)S;iUx). It 

r;— ±1 ly^fi 

turned out that the artificial configuration has a 
memory of Q = 1 for up to 10 Monte Carlo hits. 
When cooling equilibrium configurations from the 
/3 = 0.9 ensemble it never happened that they 
obtained an integer-valued topological charge, as 
depicted in Fig. |6[ 

5. Summary 

• Exact zero-modes of the overlap-Dirac operator 
are observed in the confined phase of compact 
U(l) gauge theory. 

• Agreement of the distribution of the small- 
est non-zero eigenvalue with predictions from 
random-matrix theory for fixed zero-mode degen- 
eracy is found. 

• No obvious correlations between toron charge, 
number of monopoles, number of Dirac sheets and 
the number of exact zero-modes is seen. 

• Exact zero-modes survive in the monopole part 
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Figure 6. Cooling history of typical equilibrium 
configurations at /3 = 0.9. No integer topological 
charge is found. 



of decomposed gauge field configurations, but not 
in the photon part. 

• No equilibrium configurations with integer 
topological charge are found. Cooling does not 
lead to integer values for compact U(l) gauge the- 
ory 

• Only artificially constructed configurations have 
integer topological charge. Only for ^ 10 Monte 
Carlo updates can this charge be retained after 
cooling. 

• Local spatial correlations between monopoles, 
topological charge density, and eigenfunctions of 
overlap-Dirac operator are analyzed currently to 
shed light on the confining mechanism of compact 
QED. 
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